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Properties of the ground and several collective excited states of the light nuclei '^'''''^'^'^ Mg are 
described in the framework of the angular momentum projected Generator Coordinate Method using 
the quadrupole moment as collective coordinate and the Gogny force as the effective interaction. 
The calculated excitation energies and B{E2) transition probabilities agree reasonably well with 
experiment. The results clearly indicate that both the restoration of the rotational symmetry and 
the quadrupole dynamics are key ingredients for the description of the properties of the above 
mentioned nuclei. 
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I. INTRODUCTION 



Nowadays, the region of neutron-rich nuclei around iV = 20 is the subject of active research both in the experimental 
and theoretical side. The reason is the strong experimental evidence towards the existence of quadrupole deformed 
ground states in this region. The existence of deformed ground states implies that iV = 20 is not a magic number 
for the nuclei considered, opening up the possibility for a better understanding of the mechanisms behind the shell 
structure in atomic nuclei. In addition, the extra binding energy coming from deformation can help to extend thereby 
the neutron drip line in this region far beyond what could be expected from spherical ground states. Among the 
variety of available experimental data, the most convincing evidence for a deformed ground state is found in the 
"^"^Mg nucleus where both the excitation energy of the lowest lying 2+ state |Q] and the B{E2,0^ — ^ 2+) transition 
probability ||^ have been measured. Both quantities are fairly compatible with the expectations for a rotational state. 
Theoretically, from a shell model point of view, the deformed ground states are a consequence of the lower energies 
of some intruder 2p — 2h neutron excitations into the fp shell as compared to the pure sd configuration ||^. In terms 
of the mean field picture of the nucleus, a quadrupole deformed ground state only appears after taking into account 
the zero point rotational energy correction to the mean field energy |^|-^. 

In a previous paper we have computed angular momentum projected (AMP) energy landscapes, as a function of 
the mass quadrupole moment, for the nuclei ■^^~^'^Mg and ^^~^^Si. We have found that the projection substantially 
changes the conclusions extracted from a pure mean field calculation. In all the nuclei considered, exception made of 
^'^Mg, two coexistent configurations (prolate and oblate) have been found with comparable energy indicating thereby 
that configuration mixing of states with different quadrupole intrinsic deformation had to be considered. The purpose 
of this paper is to study the effect of such configuration mixing for the nuclei •^°~^'^Mg. The Si isotopes have been 
disregarded in this work as there are indications that triaxiality effects could be relevant for the description of 
their ground states and, for the moment, our calculations are restricted to axially symmetric (K — 0) configurations. 
In our calculations we have used the Gogny force ||l^ (with the DIS parameterization ||ll|]) which is known to 
provide reasonable results for many nuclear properties like ground state deformations, moments of inertia, fission 
barrier parameters, etc, all over the periodic table. As the results presented in this paper will show, this force is 
also suited for the description of quadrupole collectivity in iV « 20 nuclei. Additional results for ^^Mg with older 
parameterizations of the Gogny force are also discussed. Finally, let us mention that similar calculations to the ones 
discussed here using the Skyrme interaction have recently been reported 



II. THEORETICAL FRAMEWORK 



To compute the properties of the ground and several collective excited states of the nuclei considered in this paper we 
have used the angular momentum projected Generator Coordinate Method (AMP-GCM) with the mass quadrupole 
moment as generating coordinate. To this end, we have used the following ansatz for the K = wave functions of 
the system 

1*^) = / dq2ofUq2o)P^o 1^(920)) (1) 
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In this expression |<p((72o)) is the set of axially symmetric (i.e. K = 0) Hartree-Fock-BogoHubov (HFB) wave functions 
generated by constraining the mass quadrupole moment to the desired values (720 = {'p{Q2o)\ — l/2(a:^ + l¥'((Z2o)) 
(please, notice that this definition is a factor 1/2 smaller than the usual definition of the intrinsic quadrupole moment). 
The intrinsic wave functions Yp{(l2^)) have been expanded in a Harmonic Oscillator (HO) basis containing 10 major 
shells and with equal oscillator lengths to make the basis closed under rotations |Q. The rotation operator in the 
HO basis has been computed using the formulas of ||l^ . 
The operator 

^00 = / cin4o(/3)e-^-'"^e-''-^«e-'^^^ (2) 

is the usual angular momentum projector with the K = Q restriction ||l^ and /i(q2o) are the "collective wave functions" 
solution of the Hill- Wheeler (HW) equation 

y dq2o'>^' (120,^20) fi{l2o) Ei J d'720-^^feo,92o)/i(92o)- (3) 

In the equation above we have introduced the projected norm Af^ {q2o, 920) — (¥'(920)1 ^oo lv(92o))' the projected 
hamiltonian kernel 7i^(920,92o) = {^{Q2o)\HPqq \(p{q2o))- As the generating states PJq |<^((72o)) are not orthogonal, 
the "collective amplitudes" fl{q2o) cannot be easily interpreted. This drawback can be easily overcome by introducing 
I p^ the so-called "natural" states 

1^0 = inlr'^' I dq2,ui{q2,)P^o l^feo)) 



which are defined in terms of the eigenstates M^(g2o) and eigenvalues of the projected norm, i.e. 
J dq2QAf^ {q2o, q2o)'^i{q2o) = "fc"fc('J'2o)- The correlated wave functions |$^) are written in terms of the natural 
states as 

k 

where the new amplitudes g'^'^ have been introduced. In terms of the amplitudes g'^'^ the collective wave functions 

giiq2o) = ^9k^ui{q2o) (4) 

k 

are defined. They are orthogonal and therefore their module squared has the meaning of a probability. The intro- 
duction of the natural states also reveals a particularity of the HW equation: if the norm has eigenvalues with zero 
value they have to be removed for a proper definition of the natural states (i.e. Hnearly dependent states are removed 
from the basis). In practical cases, in addition to the zero value eigenvalues also the eigenvalues smaller than a given 
threshold have to be removed to ensure the numerical stability of the solutions of the HW equation. In order to 
account for the fact that the mean value of the number of particles operator (<i>^ | Nj (t — tt, v) usually differs 
from the nucleus' proton and neutron numbers we have followed the usual recipe p7| , p^ of replacing the hamiltonian 
hy H — XjriNTr — Z) — — N) where Att and are chemical potentials for protons and neutrons respectively. 

Concerning the density dependent part of the Gogny force we have used the usual prescription already discussed 
in Refs. p^,p|,p^ . It amounts to use the density 

,^ ^ (y(g20)|pe-'^^" |y(g20)) /g^ 

{^{q2o)\e-^^'^W{q2,)) 

in the density dependent part of the interaction when the evaluation of {^p{q2o)\ ife^*'^'^« Wiqw)) is required in the 
calculation of the projected hamiltonian kernels. 

It has to be kept in mind that the solution of the HW equation for each value of the angular momentum / determines 
not only the ground state (cr = 1), which corresponds to the Yrast band, but also excited states ((7 = 2,3,.. .) that, 
with the set of generating wave functions used in these calculations, could correspond to solutions with a different 
deformation from the one of the ground state and/or to quadrupole vibrational bands. 

Finally, let us mention that, as the intrinsic wave functions \'^{q2o)) are determined before the projection onto 
angular momentum, the procedure described above is of the "projection after variation" (PAV) type. It is well known 



2 



that the PAV method yields the wrong moments of inertia, at least in the translational case, and a way to cure this 
deficiency is to consider a "projection before variation" (PBV) which is much more difficult to implement because the 
intrinsic wave functions have to be determined for each value of the angular momentum / using the Ritz variational 
principle on the projected energy (see for the application of PBV with small configuration spaces). To illustrate 
the consequences of the PBV method it is convenient to consider a strongly deformed intrinsic configuration \ip{q2o)) 
as in this case it is possible to obtain jl6| an approximate expression for the (PAV) projected energy Epav{I) =< 

H > — ^ ^ 2Jy^^'' where Jy is the Yoccoz (Y) moment of inertia. In this expression we recognize the rotational 

energy correction and the usual rotor- like expression for the energy of the band ^ 2Jy"^^ • shown in 

1^ (see also IH]) that starting from the projected energy and making an approximate projection before variation 
(PBV) one obtain for the energy of the rotational band the following expression EpBv{I) —< H > — ^2j^ + ^' 2^jtv^'' 
where Jty is the Thouless-Valatin (TV) moment of inertia. This implies that for the determination of the zero point 
rotational energy correction (which is very important as it can dramatically change the energy landscape as a function 
of the quadrupole moment) one has to use the Yoccoz moment of inertia (i.e. PAV is good) but for the moment of 
inertia of the band one has to use the Thouless-Valatin expression or carry out a full PBV calculation. 

Taking into account that, in the limit of strong deformation the PBV for the restoration of the rotational symmetry 
yields to the well known Self Consistent Cranking (SCC) method, a possible way to improve the AMP-GCM would be 
to consider for the intrinsic states a set of wave functions {q2o)) solution of the SCC-HFB equations for each spin 
/. However, this would lead to a triaxial projection which is extremely time consuming and also to the issue of how 
to handle configurations with (720 values close to sphericity where the SCC-HFB is no longer a good approximation 
to the PBV theory. 

In order to explore the effect of the PBV in our calculations we will restraint ourselves to perform SCC-HFB calcu- 
lations for selected configurations and compare the results with those of an AMP calculation on those configurations 
in order to extract the SCC and Yoccoz moments of inertia. The result of the comparison is that the AMP gamma 
ray energies are typically a factor 1.4 bigger than the selfconsistent ones and therefore a way to incorporate the effects 
of PBV would be to quench the bands generated by the AMP-GCM by a factor 1/1.4 « 0.7. From a physical point of 
view it is rather simple to understand why the AMP rotational band energies are higher than the SCC ones. For the 
sake of simpHcity we will concentrate on the 0+ and 2+ states. The effect of the PBV on the 0"*" state is to incorporate 
into the corresponding intrinsic state admixtures of two, four, etc quasiparticle configurations coupled to = 0. For 
the 2+ state we can also mix K = 1 and K = 2 multiquasiparticle configurations that make the variational space 
bigger and therefore leads to a higher energy gain for the 2+ state as compared to the energy gain of the 0+ state 
reducing thereby the corresponding 2+ gamma ray energy. 



III. DISCUSSION OF THE RESULTS 



A. Mean field and angular momentum projected energies 

In figure 1 we have plotted the / — Oh, 2h, 4?i, 6% and 1 — 8% projected energies E\q2o) = 'H^{q20, 920 )/A/''^ ((720, 920) 
as a function of (720 for the nuclei ^'^•^^'^'^Mg. The HFB energies have also been plotted for comparison. The projected 
energy curves can be regarded as the potential energies felt by the quadrupole collective motion and therefore give us 
indications of where the collective wave functions will be concentrated. 

Before commenting the physical contents of the curves we have to mention that, except for the I = Oh curves, 
several values around 520 = are omitted. They correspond to intrinsic configurations with a very small value of 
the norm N'^{q2Q, 920), that is, to configurations whose / — 2h, Ah, . . . contents are very small. As a consequence, the 
evaluation of the projected energies in these cases is vulnerable to strong numerical inaccuracies. Fortunately, the 
smallness of their projected norms guarantees that these configurations do not play a role in the configuration mixing 
calculation (the associated norm eigenvalues are very small) and therefore can be safely omitted. 

Coming back to the projected energy surfaces, we observe that for I = Oh and 2h a prolate and an oblate minima 
appear with almost the same energy for the nucleus ■^^Mg whereas the prolate minimum becomes deeper than the 
oblate one for 32.34^^^ Yor increasing spins either the prolate minimum becomes significantly deeper than the oblate 
one or the oblate minimum is washed out. The prolate minima are located, for all nuclei and spin values, around 
920 — that corresponds to a /3 deformation parameter of 0.4. On the other hand, the HFB energy curves show a 
behavior rather different from the I = Oh projected curves showing a spherical minimum for ■^^'^^Alg and a prolate 
one for ^^Mg. 

[FIG. 1 about here.] 
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To disentangle the relevant configurations of the intrinsic wave functions we have computed their spherical orbit 
occupancies which are given by 

i^{nlj) = {(p{q2o)\Y^ cli^mCnijm Mq2a)) (6) 

m 

where Cnijm are the annihilation operators corresponding to spherical harmonic oscillator wave functions. In the 
nucleus ^^Mg the neutron 1^(1/7/2) occupancy is zero for 520 = whereas it is almost 2 at the minimum of the 
projected energy (i.e. 920 = 1^)- The conclusion is clear, the zero point energy associated to the restoration of the 
rotational symmetry favors the configuration in which a couple of neutrons have been promoted from the sd shell to 
the /7/2 orbit. This is in good agreement with the Shell Model picture of deformation in these nuclei |Q. 

Through exhaustive mean field studies of the nucleus ^^Mg with several parameterizations of the Skyrme interaction 
1^ it has become clear that the occurrence of deformation in this nucleus is correlated to the relative position between 
the /7/2 and ^3/2 neutron orbitals. In our case (DIS parameterization of the Gogny interaction) the so-called sd — pf 
spherical shell gap for neutrons in the nucleus '^^Mg, which is given by Acf^^^-d^^^ = e/^^^ — e^gyj (with e being the 
single particle energies of the spherical configuration), takes the value 5AMeV. This value is compatible with the 
results of |^ and also with the value given in P3| . Furthermore, the /7/2— P3/2 spherical energy gap is only l.SMeV and 
therefore we expect strong quadrupole correlations between these two orbits. The values for other parameterizations 
of the Gogny force will be discussed in the last subsection. Finally, let us mention that the quantity Ae/^^^-da/a 
not well defined for the '^^Mg and ^"^Mg nuclei as in these two cases we have appreciable neutron pairing correlations 
and only the quasiparticle energies are meaningful. 



B. Angular momentum projected Generator Coordinate calculations. 

In figure 2 the collective wave functions squared |g^(g2o)|^ (see Eq. ^ for the two lowest solutions a = 1 and 
2 obtained in the AMP-GCM calculations are depicted. We also show in each panel the projected energy for the 
corresponding spin. We observe that the Of ground state wave functions of the ^^Mg and ^'^Mg nuclei contain 
significant admixtures of the prolate and oblate configurations whereas for ^'^Mg the wave function is almost completely 
located inside the prolate well. At higher spins, however, the ground state wave functions are located inside the prolate 
well in all the nuclei studied. Concerning the first excited states (cr — 2) we notice that in the nucleus ^'^Mg and 
for spins higher than zero the collective wave functions show a behavior reminiscent of a /? vibrational band: they 
are located inside the prolate wells and have a node at a 920 value near the point where the ground state collective 
wave functions attain their maximum values. Contrary to the case of a pure /3 band, the collective wave functions of 
fig. 2 are not symmetric around the node and therefore can not be considered as harmonic vibrations. On the other 
hand, the state of ^'^Mg is an admixture of prolate and oblate configurations and can not be considered as a /? 
vibrational state. The same pattern is also seen in the other two nuclei but with sHght differences: the (3 like bands 
appear at spins 4 and 6 for ^'^Mg and '^^Mg respectively. 

[FIG. 2 about here.] 

It is also worth pointing out that from the position of the tails of the collective wave functions relative to the 
projected energies (see figure caption) we can read the energy gain due to considering the quadrupole fiuctuations. The 
energy gain is maximal at / = 0/i (0.9, 1 and 0.7 MeV for ^'^Mg, ^'^Mg and ^'^Mg respectively) and quickly decreases 
with spin refiecting the narrowing of the projected wells with spin. The S{2n) separation energies are now 7.8MeV 
and 6.13MeV for ■^^Mg and ■^'^Mg respectively to be compared to the values obtained with the angular momentum 
projection |^ alone {7.65MeV and 6.39MeV) and with the experimental values of 8.056AfeT^ and 6.896MeV. 

[TABLE 1 about here.] 

In order to understand in a more quantitative way the collective wave functions just discussed it is convenient to 
analyze the quantities 

iQ2o)i ^ J dq2o\giiq2o)\^ q20, (7) 
that gives us a measure of the average deformation of the underlying intrinsic states, and 

(920)^ = / dq20 |.9^(g2o)|%20 - ((920)^) (8) 
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that serves as an estimation of the wave functions' spreading. The values of (920)0- ^^^^ = ( (^io) ) correspond 



ing to the collective wave functions of figure 2 are given in table ^ We observe that the 0^ and 23 states of ^^Mg are 
spherical (but with strong fluctuations in the (?20 degree of freedom) whereas the l\ state is deformed (/3 = 0.25). On 
the other hand, the 0^ states of ^"^Mg and ^'^Mg are deformed with (5 values of 0.16 and 0.3 respectively and have a 
S{ value rather high, possibly due to the small oblate hump. For spins higher than / = 0?i in ^^'^"^Mg and I = Ah 
in ■^°Mg the ground state band is strongly deformed. The spreading of the wave functions gets smaller for increasing 
spins as expected. The excited bands also get more deformed for increasing spin, but their f3 values never coincide 
with that of the ground state band. Obviously, their spreadings are bigger than for the ground state band. 

[TABLE 2 about here.] 

A more precise definition of the quadrupole moment for protons for each of the AMP-GCM states can be obtained 
from the results of the exact spectroscopic quadrupole moments Qa (/) for protons (no effective charge has been used) . 
The values obtained for each of the wave functions 1$^) are given in table | for the three nuclei studied and a — 1 
and 2. All the spectroscopic moments are negative indicating prolate intrinsic deformations. 

We can also compute the total intrinsic quadrupole moments from the spectroscopic ones through the formula 
(920*) (T ~ ~^^T^Qo-W^ where the K = restriction has been taken into account and also the fact that our (720 
values are, by definition, a factor 0.5 smaller than Qo- The factor A/Z is used to take into account the fact that the 
spectroscopic quadrupole moments are given in term of the proton mass distribution whereas the intrinsic quadrupole 
moments are the total ones. As can be readily observed from tabled the intrinsic quadrupole moments obtained from 
the spectroscopic ones agree rather well with the corresponding average (^20)^ 1°'* spins and deviate up to a 20% 
for spin 8h. 

[TABLE 3 about here.] 



In table [11 the energy splittings between different states and the E2 transition probabilities among them are compared 
with the available experimental data. Concerning the B{E2,0f — > 2^) transition probabilities we find a very good 
agreement with the only known experimental value and with the theoretical predictions of Utsumo et al. [Q using 
the Monte Carlo Shell Model (MCSM). The 2^ excitation energies rather nicely follow the isotopic trend but they are 
larger than the experimental values by a factor of roughly 1.5. This discrepancy could be the result of using angular 
momentum projection after variation (PAV) instead of the more complete projection before variation (PBV) that will 
require for each value of the angular momentum the calculation of the generating states from the variational principle 
on the projected energy (see section 2). Usually, the PBV method yields to rotational bands with moments of inertia 
larger than the PAV ones [p7|,po|. 

A full PBV is, unfortunately, extremely costly to implement with large configuration spaces. Therefore, to estimate 
the effect of PBV in our results, we have resorted to the selfconsistent cranking method which is an approximation to 
PBV in the limit of large deformations. We have chosen the intrinsic state with (720 = as the most representative 
configuration (it approximately corresponds to the prolate minima in all the nuclei considered) and computed the 
projected energies. In addition, selfconsistent cranking calculations with the constraints 920 — 1^ in the quadrupole 
moment and (Jx) = \/ I{I + 1) in the angular momentum have been performed. The cranking results for the exci- 
tation energies of the 2+ state are 0.548, 0.591 and 0.571 MeV for ^^Mg, ^'^Mg and ^^Mg respectively whereas the 
corresponding projected quantities are 0.753, 0.873 and 0.895 MeV. The cranking excitation energies of the 2+ state 
are a factor 0.7 smaller than the projected ones and therefore, the effect of PBV is to increase the moment of inertia 
as compared to the PAV method. If we consider the reduction factor as significative (the 520 value chosen roughly 
corresponds to the position of the maxima of the collective wave functions) and apply it to our GCM results for the 
Oi — 2^ energy differences we obtain the values 0.71, 1.02 and 1.50 MeV for ^'^Mg, ^^Mg and ^^Mg respectively. 
The new energy differences are in much better agreement with the experimental values and the MCSM results than 
the uncorrected ones. Also the corrected energy obtained for the 4^ state of '^'^Mg is in good agreement with the 
excitation energy of 2.3 MeV of a state of this nucleus which is a firm candidate to be the 4+ state belonging to the 
Yrast "rotational band" [Q. 

Although the previous estimation can be criticized in many ways we think it may serve as an indication that a full 
PBV will improve the results obtained here. Concerning the B{E2) transition probabilities, the main effect of the 
PBV will be to shift down the / — 2h, . . . projected energy curves keeping its shape mostly unaffected. Therefore, we 
do not expect big changes both in the collective wave functions 5^(^20) and in the B{E) transition probabilities that 
depend on them. 

[FIG. 3 about here.] 
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Finally, the band energy diagrams for the three nuclei considered are shown in figure 3 for states with excitation 
energies smaller than 10 MeV. For each nuclei, the bands labeled (a) and (b) correspond to the AMP-GCM result for 
the Yrast and excited bands, the band labeled (c) accounts for the experimental data in ^^Mg and ^'^Mg and for the 
MCSM result in ^'^Mg and finally, bands (d) and (e) stand for the GCM bands quenched by the factor 0.7 previously 
discussed. 



C. Results for other parameterizations of the Gogny force. 

The occurrence of quadrupole deformation in atomic nuclei is the result of the competition between two effects; 
namely, the surface energy which prevents deformation and the quantal shell effects which, depending on the nucleus, 
favor quadrupole deformation. It is therefore highly interesting to analyze the effect of these two aspects in the results 
we have obtained for the nucleus ^^Mg. To this end we have carried out calculations with two old parameterizations 
of the Gogny force; namely, the Dl and Dl' parameterizations The Dl parameterization was the one originally 
proposed by Gogny and the only difference with Dl' is the spin-orbit strength which is smaller for Dl. As a result 
one can expect that Dl will lead to a higher Aef^^^-d3/2 = ^hii ~ ^^3/2 energy gap than Dl' as it turns out to be the 
case: the value of Ae/^^2-d3/2 6.37MeT^ for Dl and 5.37MeV for Dl'. On the contrary, the value of Aef^^^-p^^^for 
Dl gets reduced from the 1.91MeV we obtain for Dl' to the value 1.56MeV. On the other hand, DIS has the same 
spin-orbit strength as Dl' (the values of Ae/^^^-da/a and Acf^^^-p^^^ given in the previous paragraph for Dl' are very 
close to those of DIS given in a previous subsection) but its surface energy coefficient is smaller than in Dl'. The 
need for a reduction of the surface energy coefficient in Dl' was evident when the fission barriers for ^'^^Pu ||ll) were 
computed with the Gogny force: they came out too high and the new DIS parameterization was proposed to cure 
this deficiency of the former Dl' parameterization. 

In figure 4 we have plotted the HFB energy curves (left panel) and the AMP energies for / = (middle panel) and 
1 — 2 (right panel) for the three parameterizations of the Gogny force just mentioned. We observe that the results 
obtained for DIS and Dl' are, apart from the overall AMeV shift, very similar. This similarity is a clear indication that 
the value of the surface energy parameter has no infiuence on the results. The HFB result for Dl shows a shoulder at 
120 = which is located much higher in energy than the corresponding shoulder for DIS and Dl'. As a consequence, 
the 7 = projected energy curve obtained with Dl shows a very shallow minimum at (720 = 0.56. However, the 
1 = 2 projected energy curves are very similar for the three parameterizations. The differences found between the Dl 
results and the ones with the two other parameterizations clearly indicate the sensitivity of the quadrupole properties 
of ■^^Mg to the relative position of the orbits involved. 

[FIG. 4 about here.] 

Finally, we have carried out the AMP-GCM calculation for the Dl and Dl' parameterizations of the force and the 



most important quantities obtained are summarized in table IV. As expected from the projected energy curves of figure 



4 we obtain a rather small average quadrupole moment (920)0- ^'^^ cr = 1 and I = Oh with the Dl parameterization 
and bigger ones for the two other parameterizations. However, the (920)^ for f = 1 and 1 = 2% are rather similar in 
the three cases. The smaller value of (^20)1 the Dl parameters gets refiected in a much smaller B{E2) transition 
probabilities than for the two other parameterizations. Finally, the excitation energy of the 2f state with respect to 
the ground state turns out to be significantly bigger for Dl than for the other parameter sets, being the results of Dl' 
and DIS in reasonable agreement. The final conclusion of this comparison is that the energy gap Ae/^^2-rf3/2 seems 
to be a relevant parameter in order to reproduce the properties of '^^Mg. 

[TABLE 4 about here.] 



IV. CONCLUSIONS 



In conclusion, we have performed angular momentum projected Generator Coordinate Method calculations with 
the Gogny interaction DIS and the mass quadrupole moment as generating coordinate in order to describe rotational 
like states in the nuclei ^^Mg, ^^Mg and ^"^Mg. We obtain a very well deformed ground state in ^'^Mg, a fairly 
deformed ground state in ^^Mg and a spherical ground state in ■^^Mg. In the three nuclei, states with spins higher or 
equal I = Afi are deformed. The intraband B{E2) transition probabilities agree well with the available experimental 
data and results from shell model like calculations. The 2+ excitation energies follow the isotopic trend but come out 
a factor 1.5 too high as compared with the experiment. We attribute the discrepancy to the well known deficiency of 
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Projection After Variation calculations of providing small moments of inertia. However, we consider the agreement 
with experiment to be remarkable taking into account that the same force used in this calculation is also able to give 
reasonable values for such different quantities as fission barrier heights, moments of inertia of superdeformed bands, 
the energy of octupole vibrations, etc in heavy nuclei. The sensitivity of the results to other parameterization of the 
Gogny interaction is also analyzed and the conclusion is that the Dl parameter set fails to reproduce the properties 
of ^^Mg being the spin-orbit strength the responsible for such failure. 
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FIG. 1. The HFB (dashed line) and angular momentum projected energies (7 = 0?i,2fi, . . . ,8fi) for the nuclei considered. 
See text for further comments. 
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q2o (b) q2o (b) q2o (b) 

FIG. 2. The collective amplitudes \gi{q2o)\'' (thick lines) for a = 1 (full) and 2 (dashed) and spin values of / = OR, . . . , S/i 
for the nuclei ^"Mg, '^^Mg and ^'^Mg. The projected energy curve for each spin is also plotted (thin line). The y-axis scales 
are in energy units and always span an energy interval of 13 MeV (minor ticks are 0.5 MeV apart). The collective wave 
functions \gi{q2o)f have also been plotted against the energy scale after a proper scaling and shifting, that is, the quantity 
Ei + 15 X |.9^(g2o)P is the one actually plotted. With this choice of scales we can read from the figure the energy gain due to 
the quadrupole fluctuations by considering the position of the wave functions' tail relative to the projected curve. 
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FIG. 3. Collective bands for the three nuclei studied. Bands a) and b) correspond to the AMP-GCM results for the ground 
and first excited band, c) is the experimental band (in the case of ^^Mg the MCSM prediction of Q has been used). Finally, 
bands d) and e) are the AMP-GCM results quenched by the factor 0.7 discussed in the text 
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^^Mg 


i 




^1 








^1 

















0.091 


0.558 


0.626 


0.685 


0.436 


0.692 


0.396 


0.601 


0.788 


0.691 


0.440 


0.723 


2 


0.579 


0.588 


0.092 


0.750 


0.885 


0.482 


0.393 


0.859 


1.052 


0.455 


0.644 


0.658 


4 


0.962 


0.387 


0.215 


0.716 


1.012 


0.388 


1.041 


0.723 


1.136 


0.387 


0.819 


0.573 


6 


1.087 


0.300 


0.581 


0.557 


1.084 


0.363 


1.264 


0.554 


1.188 


0.354 


0.860 


0.528 


8 


1.131 


0.289 


1.470 


0.562 


1.151 


0.368 


1.293 


0.515 


1.226 


0.332 


0.926 


0.560 



TABLE 1. The average intrinsic quadrupole moment (^20)^ ^^'^ fluctuations = y (920)^ barns for the three nuclei 
considered. 
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a = 1 


(7 = 2 




2 


4 


6 


8 


2 


4 


6 


8 




-13.79 


-27.01 


-32.43 


-35.36 


-3.11 


-10.07 


-21.48 


-38.35 




-19.15 


-26.31 


-30.09 


-31.75 


-8.63 


-23.07 


-27.22 


-29.01 




-20.78 


-27.59 


-31.27 


-33.70 


-15.16 


-21.58 


-25.34 


-26.1 



TABLE II. Spectroscopic quadrupole moments in efn? for 7 = 2?i, 47i, 67i and 87i and tr = 1 and 2 for the three nuclei 
considered in this paper. 
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Calc. Energies (MeV) 


Exp. 


Calc. B{E2)e^fm'^ 


Exp. 




0+-2+ 




2+ - 2+ 


0+-2+ 


0+-2+ 


0+^2+ 


0+^2+ 


0+-2+ 




2.15 


2.30 


1.60 


1.482 


229 


3 


218 


300(*) 




1.46 


1.77 


3.35 


0.885 


395 


3.4 


199 


454±78 




1.02 


2.35 


3.31 


0.75(*) 


525 





290 


580(*) 



TABLE III. Calculated and experimental results for excitation energies and B{E2,0'^^ — > 2'^^) transition probabilities. In 
the experimental data columns values marked with an (*) correspond to Monte Carlo Shell Model results taken from Ref. Ipif . 
The experimental data for the excitation energies have been taken from Q for the ^^Mg nucleus and from j25| for '^^ Mg. The 
B{E2) transition probability has been taken from jij. 
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(920)? 


(320)? 


B(i?2,0+ ^2+) 


^0+-2+ 


^6/7/2-''3/2 


^^/7/2-P3/2 


Dl 


0.185 


0.785 


138 


2.25 


6.37 


1.56 


Dl' 


0.381 


0.869 


299 


1.67 


5.37 


1.91 


DIS 


0.436 


0.885 


385 


1.46 


5.37 


1.80 



TABLE IV. Results of the AMP-GCM calculations for '^'^Mg and the parameterizations Dl, Dl' and DIS of the Gogny 
interaction. The average quadrupole moments (92o)ct ground state band and spins and 2 are given, in barns, in 

the first two columns. In the third column the B{E2) transition probabilities in e^fmf^ are given. In the fourth column the 
excitation energy of the 2+ state with respect to the ground state is given in MeV. Finally, in the last two columns the energy 
gaps A£/^^2-ci3/2 and Aef^^^-p^^^ are given in MeV. 
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